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The characterization of the dynamics of quantum systems is a task of both fundamental and practical impor- 
tance. A general class of methods which have been developed in quantum information theory to accomplish this 
task is known as quantum process tomography (QPT). In an earlier paper [M. Mohseni and D. A. Lidar, Phys. 
Rev. Lett. 97, 170501 (2006)] we presented a new algorithm for Direct Characterization of Quantum Dynamics 
(DCQD) of two-level quantum systems. Here we provide a generalization by developing a theory for direct and 
complete characterization of the dynamics of arbitrary quantum systems. In contrast to other QPT schemes, 
DCQD relies on quantum error-detection techniques and does not require any quantum state tomography. We 
demonstrate that for the full characterization of the dynamics of n d-level quantum systems (with d a power 
of a prime), the minimal number of required experimental configurations is reduced quadratically from d in in 
separable QPT schemes to d 2n in DCQD. 
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I. INTRODUCTION 

The characterization of quantum dynamical systems is a 
fundamental problem in quantum physics and quantum chem- 
istry. Its ubiquity is due to the fact that knowledge of quantum 
dynamics of (open or closed) quantum systems is indispens- 
able in prediction of experimental outcomes. In particular, 
accurate estimation of an unknown quantum dynamical pro- 
cess acting on a quantum system is a pivotal task in coherent 
control of the dynamics, especially in verifying/monitoring 
the performance of a quantum device in the presence of de- 
coherence. The procedures for characterization of quantum 
dynamical maps are traditionally known as quantum process 
tomography (QPT) GL0!]. 

In most QPT schemes the information about the quan- 
tum dynamical process is obtained indirectly. The quantum 
dynamics is first mapped onto the state(s) of an ensemble 
of probe quantum systems, and then the process is recon- 
structed via quantum state tomography of the output states. 
Quantum state tomography is itself a procedure for identify- 
ing a quantum system by measuring the expectation values 
of a set of non-commuting observables on identical copies 
of the system. There are two general types QPT schemes. 
The first is Standard Quantum Process Tomography (SQPT) 
0, IU> IH . In SQPT all quantum operations, including prepa- 
ration and (state tomography) measurements, are performed 
on the system whose dynamics is to be identified (the "prin- 
cipal" system), without the use of any ancillas. The SQPT 
scheme has already been experimentally demonstrated in a 
variety of systems including liquid-state nuclear magnetic res- 
onance (NMR) U lip, optical U [13], atomic JH, and 
solid-state systems 1121 . The second type of QPT scheme 
is known as Ancilla-Assisted Process Tomography (AAPT) 
E [Mill dl- In AAPT one makes use of an ancilla (auxil- 
iary system). First, the combined principal system and ancilla 
are prepared in a "faithful" state, with the property that all in- 
formation about the dynamics can be imprinted on the final 
state lfl3l n~5l [Hill . The relevant information is then extracted 
by performing quantum state tomography in the joint Hilbert 



space of system and ancilla. The AAPT scheme has also been 
demonstrated experimentally lfl5i [l7tl . The total number of 
experimental configurations required for measuring the quan- 
tum dynamics of n d-level quantum systems (qudits) is d An 
for both SQPT and separable AAPT, where separable refers 
to the measurements performed at the end. This number can 
in principle be reduced by utilizing non-separable measure- 
ments, e.g., a generalized measurement j^. However, the 
non-separable QPT schemes are rather impractical in physi- 
cal applications because they require many-body interactions, 
which are not experimentally available or must be simulated 
at high resource cost (H). 

Both SQPT and AAPT make use of a mapping of the 
dynamics onto a state. This raises the natural question of 
whether it is possible to avoid such a mapping and instead per- 
form a direct measurement of quantum dynamics, which does 
not require any state tomography. Moreover, it seems reason- 
able that by avoiding the indirect mapping one should be able 
to attain a reduction in resource use (e.g., the total number 
of measurements required), by eliminating redundancies. In- 
deed, there has been a growing interest in the development 
of direct methods for obtaining specific information about the 
states or dynamics of quantum systems. Examples include the 
estimation of general functions of a quantum state [18], detec- 
tion of quantum entanglement lfl9tl . measurement of nonlinear 
properties of bipartite quantum states 1I20I1 . reconstruction of 
quantum states or dynamics from incomplete measurements 
12111 . estimation of the average fidelity of a quantum gate or 
process |22[ [2311 . and universal source coding and data com- 
pression Il24ll . However, these schemes cannot be used di- 
rectly for a complete characterization of quantum dynamics. 
In Ref. I25I1 we presented such a scheme, which we called 
"Direct Characterization of Quantum Dynamics" (DCQD). 

In trying to address the problem of direct and complete 
characterization of quantum dynamics, we were inspired by 
the observation that quantum error detection (QED) |Q]] pro- 
vides a means to directly obtain partial information about the 
nature of a quantum process, without ever revealing the state 
of the system. In general, however, it is unclear if there is 
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FIG. 1: Schematic of DCQD for a single qubit, consisting of Bell- 
state preparations, application of the unknown quantum map, £, and 
Bell-state measurement (BSM). 



a fundamental relationship between QED and QPT, namely 
whether it is possible to completely characterize the quantum 
dynamics of arbitrary quantum systems using QED. And, pro- 
viding the answer is affirmative, how the physical resources 
scale with system size. Moreover, one would like to under- 
stand whether entanglement plays a fundamental role, and 
what potential applications emerge from such a theory linking 
QPT and QED. Finally, one would hope that this approach 
may lead to new ways of understanding and/or controlling 
quantum dynamical systems. We addressed these questions 
for the first time in Ref. ll25ll by developing the DCQD algo- 
rithm in the context of two-level quantum systems. In DCQD 
- see Fig. [TJ- the state space of an ancilla is utilized such that 
experimental outcomes from a Bell-state measurement pro- 
vide direct information about specific properties of the under- 
lying dynamics. A complete set of probe states is then used to 
fully characterize the unknown quantum dynamics via appli- 
cation of a single Bell-state measurement device ll3i l25ll . 

Here we generalize the theory of Ref. l25ll to arbitrary 
open quantum systems undergoing an unknown, completely- 
positive (CP) quantum dynamical map. In the generalized 
DCQD scheme, each probe qudit (with d prime) is initially 
entangled with an ancillary qudit system of the same dimen- 
sion, before being subjected to the unknown quantum pro- 
cess. To extract the relevant information, the corresponding 
measurements are devised in such a way that the final (joint) 
probability distributions of the outcomes are directly related 
to specific sets of the dynamical superoperator's elements. A 
complete set of probe states can then be utilized to fully char- 
acterize the unknown quantum dynamical map. The prepara- 
tion of the probe systems and the measurement schemes are 
based on QED techniques, however, the objective and the de- 
tails of the error-detection schemes are different from those 
appearing in the protection of quantum systems against deco- 
herence (the original context of QED). More specifically, we 
develop error-detection schemes to directly measure the co- 
herence in a quantum dynamical process, represented by off- 
diagonal elements of the corresponding superoperator. We ex- 
plicitly demonstrate that for characterizing a dynamical map 
on n qudits, the number of required experimental configura- 
tions is reduced from d in , in SQPT and separable AAPT, to 
d 2n in DCQD. A useful feature of DCQD is that it can be 
efficiently a ppli ed to partial characterization of quantum dy- 
namics I25I I26I1 . For example, it can be used for the task of 
Hamiltonian identification, and also for simultaneous deter- 
mination of the relaxation time T\ and the dephasing time T2. 

This paper is organized as follows. In Sec. [TTJ we provide a 
brief review of completely -positive quantum dynamical maps, 



and the relevant QED concepts such as stabilizer codes and 
normalizers. In Sec. [Till we demonstrate how to determine 
the quantum dynamical populations, or diagonal elements of 
a superoperator, through a single (ensemble) measurement. In 
order to further develop the DCQD algorithm and build the 
required notations, we introduce some lemmas and definitions 
in Sec. HVI and then we address the characterization of quan- 
tum dynamical coherences, or off-diagonal elements of a su- 
peroperator, in Sec. [V] In Sec. I VII we show that measurement 
outcomes obtained in Sec. [VI provide d 2 linearly independent 
equations for estimating the coherences in a process, which is 
in fact the maximum amount of information that can be ex- 
tracted in a single measurement. A complete characterization 
of the quantum dynamics, however, requires obtaining d A in- 
dependent real parameters of the superoperator (for non-trace 
preserving maps). In Sec. IVHI we demonstrate how one can 
obtain complete information by appropriately rotating the in- 
put state and repeating the above algorithm for a complete 
set of rotations. In Sec. I Villi and IIX1 we address the general 
constraints on input stabilizer codes and the minimum num- 
ber of physical qudits required for the encoding. In Sec. [X] 
and Sec. [XI] we define a standard notation for stabilizer and 
normalizer measurements and then provide an outline of the 
DCQD algorithm for the case of a single qudit. For conve- 
nience, we provide a brief summary of the entire DCQD al- 
gorithm in Sec. IXIII We conclude with an outlook in Sec- 
tion [XlII] In Appendix lAl we generalize the scheme for arbi- 
trary open quantum systems. For a discussion of the experi- 
mental feasibility of DCQD see Ref. l25ll . and for a detailed 
and comprehensive comparison of the required physical re- 
sources in different QPT schemes see Ref. ^ . 



II. PRELIMINARIES 

In this section we introduce the basic concepts and nota- 
tion from the theory of open quantum system dynamics and 
quantum error detection, required for the generalization of the 
DCQD algorithm to qudits. 



A. Quantum Dynamics 

The evolution of a quantum system (open or closed) can, 
under natural assumptions, be expressed in terms of a com- 
pletely positive quantum dynamical map £, which can be rep- 
resented as 0]] 

d 2 -i 

S (P)= J2 XmnEmpEl (1) 
m,n— 

Here p is the initial state of the system, and the {E m } are a set 
of (error) operatorbasis elements in the Hilbert-Schmidt space 
of the linear operators acting on the system. I.e., any arbitrary 
operator acting on a d-dimensional quantum system can be ex- 
panded over an orthonormal and unitary error operator basis 
{i? , Ex,..., i? (i 2„ 1 }, where E = I and tv(EjEj) = d5 tj 
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ll27ll . The {xmn} are the matrix elements of the superopera- 
tor x, or "process matrix", which encodes all the information 
about the dynamics, relative to the basis set {E m } |Q]]. For an 
?i-qudit system, the number of independent matrix elements 
in x is d 4n for a non-trace-preserving map and d 4n — d 2n for 
a trace-preserving map. The process matrix x is positive and 
Trx < 1. Thus x can t> e thought of as a density matrix in 
the Hilbert-Schmidt space, whence we often refer to its diag- 
onal and off-diagonal elements as "quantum dynamical popu- 
lation" and "quantum dynamical coherence", respectively. 

In general, any successive operation of the (error) opera- 
tor basis can be expressed as EiEj = J^k ujl, -'' k Ek, where 
i, j, k = 0, 1, . . . , d 2 — 1. However, we use the "very nice (er- 
ror) operator basis" in which EiEj = co^E^j, det Ei = 1, 
uj % -i is a dth root of unity, and the operation * induces a 
group on the indices (27 1. This provides a natural general- 
ization of the Pauli group to higher dimensions. Any ele- 
ment Ei can be generated from appropriate products of Xd 
and Z d , where X d \k) = \k + l) , Z d \k) = uj k [k), and 
X d Z d = uJ~ 1 Z d Xd l27ll28ll . Therefore, for any two elements 

Ei={a, 9 *} = " a X 9 d Z p d and E j={a ,, q ,, p , } = u, a ' X% Z p d ' 
(where < q,p < d) of the single-qudit Pauli group, we 
always have 

E i E j =uj pq '-' IP 'E j E i , (2) 

where 

pq' — qp 1 = k (mod d). (3) 

The operators Ei and Ej commute iff k — 0. Henceforth, 
all algebraic operations are performed in mod(e?) arithmetic, 
and all quantum states and operators, respectively, belong to 
and act on a <i-dimensional Hilbert space. For simplicity, from 
now on we drop the subscript d from the operators. 

B. Quantum Error Detection 

In the last decade the theory of quantum error correction 
(QEC) has been developed as a general method for detect- 
ing and correcting quantum dynamical errors acting on multi- 
qubit systems such as a quantum computer |Q]]. QEC consists 
of three steps: preparation, quantum error detection (QED) 
or syndrome measurements, and recovery. In the preparation 
step, the state of a quantum system is encoded into a subspace 
of a larger Hilbert space by entangling the principal system 
with some other quantum systems using unitary operations. 
This encoding is designed to allow detection of arbitrary er- 
rors on one (or more) physical qubits of a code by performing 
a set of QED measurements. The measurement strategy is to 
map different possible sets of errors only to orthogonal and 
undeformed subspaces of the total Hilbert space, such that the 
errors can be unambiguously discriminated. Finally the de- 
tected errors can be corrected by applying the required unitary 
operations on the physical qubits during the recovery step. A 
key observation relevant for our purposes is that by perform- 
ing QED one can actually obtain partial information about the 
dynamics of an open quantum system. 



For a qudit in a general state \<p c ) in the code space, 
and for arbitrary error basis elements E m and E n , the 
Knill-Laflamme QEC condition for degenerate codes is 
(</> c | E^Em \(f> c ) = a nm , where a nm is a Hermitian matrix 
of complex numbers Q|]. For nondegenerate codes, the QEC 
condition reduces to (</> c | E^Em \(f> c ) = 8 nm ; i.e., in this case 
the errors always take the code space to orthogonal subspaces. 
The difference between nondegenerate and degenerate codes 
is illustrated in Fig. [2] In this work, we concentrate on a 
large class of error-correcting codes known as stabilizer codes 
112911 ; however, in contrast to QEC, we restrict our attention 
almost entirely to degenerate stabilizer codes as the initial 
states. Moreover, by definition of our problem, the recov- 
ery/correction step is not needed or used in our analysis. 

A stabilizer code is a subspace He of the Hilbert space of 
n qubits that is an eigenspace of a given Abelian sub grou p S 
of the n-qubit Pauli group with the eigenvalue +1 111 l29ll . In 
other words, for \<f> c ) 6 He an d Si £ S, we have Si \ 4> c ) = 
1 C ), where Si's are the stabilizer generators and [Si, Sj] = 0. 
Consider the action of an arbitrary error operator E on the 
stabilizer code \<fi c ), E \ <f> c }- The detection of such an error 
will be possible if the error operator anticommutes with (at 
least one of) the stabilizer generators, {Si, E} = 0. I.e, by 
measuring all generators of the stabilizer and obtaining one or 
more negative eigenvalues we can determine the nature of the 
error unambiguously as: 

Si(E \cf> c }) = -E(Si \<j> c }) = -{E \4> c )). 

A stabilizer code [n,k,d c ] represents an encoding of k log- 
ical qudits into n physical qudits with code distance d c , such 
that an arbitrary error on any subset of t = (d c — l)/2 or 
fewer qudits can be detected by QED measurements. A sta- 
bilizer group with n — k generators has d n ~ k elements and 
the code space is d fe -dimensional. Note that this is valid when 
d is a power of a prime ll28ll . The unitary operators that pre- 
serve the stabilizer group by conjugation, i.e., USU^ = S, 
are called the normalizer of the stabilizer group, N(S). Since 
the normalizer elements preserve the code space they can be 
used to perform certain logical operations in the code space. 
However, they are insufficient for performing arbitrary quan- 
tum operations [U] . 

Similarly to the case of a qubit l25ll . the DCQD algorithm 
for the case of a qudit system consists of two procedures: (i) 
a single experimental configuration for characterization of the 
quantum dynamical populations, and (ii) d 2 — 1 experimental 
configurations for characterization of the quantum dynamical 
coherences. In both procedures we always use two physical 
qudits for the encoding, the principal system A and the ancilla 
B, i.e., n = 2. In procedure (i) - characterizing the diagonal 
elements of the superoperator - the stabilizer group has two 
generators. Therefore it has d 2 elements and the code space 
consists of a single quantum state (i.e., k = 0). In procedure 
(ii) - characterizing the off-diagonal elements of the superop- 
erator - the stabilizer group has a single generator, thus it has 
d elements, and the code space is two-dimensional. That is, 
we effectively encode a logical qudit (i.e., k = 1) into two 
physical qudits. In next sections, we develop the procedures 
(i) and (ii) in detail for a single qudit with d being a prime, 
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FIG. 2: A schematic diagram of Quantum Error Detection (QED). The projective measurements corresponding to eigenvalues of stabilizer 
generators are represented by arrows. For a non-degenerate QEC code, after the QED, the wavefunction of the multiqubit system collapses 
into one of the orthogonal subspaces each of which is associated with a single error operator. Therefore, all errors can be unambiguously 
discriminated. For degenerate codes, by performing QED the codespace also collapses into a set orthogonal subspaces. However, each 
subspace has multiple degeneracies among k error operators in a subset of the operator basis, i.e., {E m }^ n _ 1 C {Ei}^^ 1 . In this case, one 
cannot distinguish between different operators within a particular subset {-E m } ° =1 . 



and in the appendix [A] we address the generalization to sys- 
tems with d being an arbitrary power of a prime. 

III. CHARACTERIZATION OF QUANTUM DYNAMICAL 
POPULATION 

To characterize the diagonal elements of the superoperator, 
or the population of the unitary error basis, we use a non- 
degenerate stabilizer code. We prepare the principal qudit, 
A, and an ancilla qudit, B, in a common +1 eigenstate \<p c ) 
of the two unitary operators E A E B and E A E B , such that 
[E A Ef ,E A Ef,] = (e.g. X A X B and Z A (Z B ) d ~ 1 ). There- 
fore, simultaneous measurement of these stabilizer generators 



at the end of the dynamical process reveals arbitrary single 
qudit errors on the system A. The possible outcomes depend 
on whether a specific operator in the operator-sum represen- 
tation of the quantum dynamics commutes with E A E B and 
E A E B , with the eigenvalue +1, or with one of the eigenval- 
ues ui, lu 2 , . . . , u> x . The projection operators correspond- 
ing to outcomes u k and uj k , where k,k' = 0, 1, . . . , d — 1, 
have the form P k = u~ lk {E A Ef) 1 and P k , = 

1 Ef'=o UJ ^ l ' k ' i E v E j> f ■ The J oint probability distribution 
of the commuting Hermitian operators P k and P k > on the out- 
put state £ 0) = Em,nX™« E mPEl, where p = \<j> c ) (<f) c \ , 
is: 



d 2 -l d-ld-1 

Tr[P k P k ,£(p)] = ^ £ X mn J2zZ u ~ lku ~ l ' k '^^n(Ef) l (E^ l 'E m (Ef)\Ep) l 'p]. 

m,n=0 1=0 l'=0 

Using E l E m = u> im E m Ei and the relation (E A Ef) l (E A Ef, ) v p = p, we obtain: 

d 2 -l 

Tr[P k P k ,£(p)] = - 2 £ X mnY.zZ^ k)lujK ^ k ' )l ' 6 ^ 

m,n=0 1=0 l'=0 

I 



where we have used the QED condition for nondegenerate 
codes: 

Tr[E\E m p] = {<t> c \ E\E m \<f> c ) = S mn , 

i.e., the fact that different errors should take the code space 
to orthogonal subspaces, in order for errors to be unam- 
biguously detectable, see Fig. [3] Now, using the discrete 
Fourier transform identities X^/tTrj 1 u/ lm_fc ^ = dSi mtk and 



J2i'=o U3 ^' m k '^ 1 ' = ^i' m ,k', we obtain: 
d 2 -i 

Tr[P k P k ,£(p)} = 2J Xmrn k m ,k^i' m ,k' = Xmomo- (4) 
m=0 

Here, mo is defined through the relations i mo = k and 
i' = k', i.e., E m<1 is the unique error operator that anti- 
commutes with the stabilizer operators with a fixed pair of 
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eigenvalues uu k and oj k corresponding to the experimental 
outcomes k and k'. Since each P k and Py operator has d 
eigenvalues, we have d 2 possible outcomes, which gives us 
d 2 linearly independent equations. Therefore, we can char- 
acterize all the diagonal elements of the superoperator with a 
single ensemble measurement and 2d detectors. 

In order to investigate the properties of the pure state \<j) c ), 
we note that the code space is one-dimensional (i.e., it has 
only one vector) and can be Schmidt decomposed as \<p c ) = 
12k=o \k) A \k) B , where A& are non-negative real numbers. 
Suppose Z\k) = ui k \k); without loss of generality the two 
stabilizer generators of \<p c ) can be chosen to be {X A X B ) q 
and [Z A (Z B ) d - 1 ]P. We then have (0 C | (X A X B ) q \<j> c ) = 1 
and \<j> c \ [Z A {Z B ) d - 1 ]P \<t> c ) = 1 for any q and p, where < 
q,p < d. This results in the set of equations X)fc=o ^k^k+q = 
1 for all q, which have only one positive real solution: Ao = 
Ai = ... = Afc = 1/vb; i.e., the stabilizer state, \4> c ), is 
a maximally entangled state in the Hilbert space of the two 
qudits. 

In the remaining parts of this paper, we first develop an 
algorithm for extracting optimal information about the dy- 
namical coherence of a d-level quantum system (with d be- 
ing a prime), through a single experimental configuration, in 
Sec.[lV][V]and[Vl] Then, we further develop the algorithm to 
obtain complete information about the off-diagonal elements 
of the superoperator by repeating the same scheme for differ- 
ent input states, Sec. lVIII In Sec. [A] we address the general- 
ization of the DCQD algorithm for qudit systems with d being 
a power of a prime. In the first step, in the next section, we 
establish the required notation by introducing some lemmas 
and definitions. 



IV. BASIC LEMMAS AND DEFINITIONS 

Lemma 1 Let < q,p, q 1 ,p' < d, where d is prime. Then, 
for given q, p, q' and fc(mod d), there is a unique p' that solves 

pq' — qp' = k (mod d). 

Proof: We have pq' — qp' = k (mod d) = k + td, where 
t is an integer. The possible solutions for p' are indexed by t 
as p'(t) = (pq' — k — td) I 'q. We now show that if p'(ti) is 
a solution for a specific value t\, there exists no other integer 
t-z 7^ ii such that p'fa) is another independent solution to this 
equation, i.e., p'fe) ^ p'(£i)(mod d). First, note that if p'fa) 
is another solution then we havep'(£i) = p'(£2)+(£2 — ti)d/q. 
Since d is prime, there are two possibilities: a) q divides (t? — 
t\), then (£2 — t\)d/q = ±nd, where n is a positive integer; 
therefore we have p'fa) = f/(£i)(mod d), which contradicts 
our assumption that ^'(£2) is an independent solution from 



p'(t\). b) q does not divide (£2 — £1), then (£2 — ti)d/q is not 
a integer, which is unacceptable. Thus, we have £2 = £1, i.e., 
the solution p'(t) is unique. 

Note that the above argument does not hold if d is not 
prime, and therefore, for some q' there could be more than 
one p' that satisfies pq' — qp' = k (mod d). In general, the 
validity of this lemma relies on the fact that Z<j is a field only 
for prime d. 

Lemma 2 For any unitary error operator basis Ei acting on 
a Hilbert space of dimension d, where d is a prime and 
i = 0, 1, . . . , d 2 — 1, there are d unitary error operator basis 
elements, Ej, that anticommute with Ei with a specific eigen- 
value u k , i.e., EiEj = uj k EjEi, where k = 0, . . . , d — 1. 

Proof: We have EiEj = uo pq '- qp ' EjE u where < 
q,p,q',p' < d, and pq' — qp' = k (mod d). Therefore, for 
fixed q, p, and k (mod d) we need to show that there are d 
solutions (q',p'). According to LemmaQ] for any q' there is 
only one p' that satisfies pq' — qp' = k (mod d) ; but g'can 
have d possible values, therefore there are d possible pairs of 

(«V). 

Definition 1 We introduce d different subsets, WL k — 
0, 1, . . . , d — 1, of a unitary error operator basis {Ej} (i.e. 
W'l C {Ej}). Each subset contains d members which all 
anticommute with a particular basis element Ei , where i = 
0,1, ... ,d 2 — 1, with fixed eigenvalue cu k . The subset Wq 
which includes Eq and Ei is in fact an Abelian subgroup of 
the single-qudit Pauli group, G\. 



V. CHARACTERIZATION OF QUANTUM DYNAMICAL 
COHERENCE 

For characterization of the coherence in a quantum dynam- 
ical process acting on a qudit system, we prepare a two-qudit 
quantum system in a non-separable eigenstate \4>ij) of a uni- 
tary operator SV, = E A E B . We then subject the qudit A to 
the unknown dynamical map, and measure the sole stabilizer 
operator Sij at the output state. Here, the state \<f>ij) is in 
fact a degenerate code space, since all the operators E A that 
anticommute with E A , with a particular eigenvalue ui k , per- 
form the same transformation on the code space and cannot 
be distinguished by the stabilizer measurement. If we express 
the spectral decomposition of Sij = E A E B as Sij = J^k 
ui k Pk, the projection operator corresponding to the outcome 
Lu k can be written as P k = ~ YaZq uj- lk (E A E B Thepost- 
measurement state of the system, up a normalization factor, 
will be: 



p k £( P )Pk = 1 ]T fali^-V'^fj'v!^^)'']- 

m,n=0 1=0 i'=0 



FIG. 3: A diagram of the error-detection measurement for estimating quantum dynamical population. The arrows represent the projection 
operators P^Py corresponding to different eigenvalues of the two stabilizer generators S and S'. These projective measurements result in 
a projection of the wavefunction of the two-qudit systems, after experiencing the dynamical map, into one of the orthogonal subspaces each 
of which is associated to a specific error operator basis. By calculating the joint probability distribution of all possible outcomes, P k Py, for 
k, k' = 0, . . . , d, we obtain all d 2 diagonal elements of the superoperator in a single ensemble measurement. 

Using the relations E t E m = u im E m Ei, E\e\ = uj-^e\eI and {E^Ef) l p{E^Ef ] f = p we have: 

1 d-l d-l d 2 -l 

P k E(p)P k = -pYt^'^Y,"^ 1 ' E XmnEmpEl 
1=0 l'=0 m,n=0 

Now, using the discrete Fourier transform properties ^2f = Q a/ 4m ~ fc ^ = G?<5; mi fc and X)f'=o uj^ 1 ^ 1 = dSi n: k , we obtain: 

P k £ (p)P k = Xmm E*pE$ + £ ( Xm n E*pE? + X * mn E£pE$ ) . (5) 
m m<n 

i 



Here, the summation runs over all E^ and E^ that belong 
to the same W\; see Lemma [2] I.e., the summation is over 
all unitary operator basis elements E^ and Ef that anti- 
commute with Ef with a particular eigenvalue ui k . Since 
the number of elements in each Wk is d, the state of the 
two-qudit system after the projective measurement comprises 
d + 2[d(d — l)/2] = d 2 terms. The probability of getting the 
outcome uj k is: 

TT[P k £(p)} = Xmm + 2 Re iXrnn Tr(E^ E£ P )}. 

(6) 

Therefore, the normalized post-measurement states are pk = 
Pk£ (p)Pk/Tr[Pk£ (p)]- These d equations provide us with in- 
formation about off-diagonal elements of the superoperator iff 



Tr[(E^Y E^p] 7^ 0. Later we will derive some general prop- 
erties of the state p such that this condition can be satisfied. 

Next we measure the expectation value of any other uni- 
tary operator basis element T rs = E^Ef on the output state, 
such that Ef ^ I, Ef ^ I, T rs e N(S) and T rs ^ {S tj ) a , 
where < a < d. Let us write the spectral decomposition 

of T rs as T rs = '^2 l uj k P^. The joint probability distribution 

fc' 

of the commuting Hermitian operators Pk and Pk> on the out- 
put state £(p) is Ti[Pk> Pk£(p)}- The average of these joint 
probability distributions of Pk and Py over different values 
of fc' becomes: w fe 'Tr[P fc ,P fe £Go)] = Ti[T rs P k £{p)\ = 
Tr (T rs pk), which can be explicitly written as: 



Tv(T rsPk ) = £ Xmm Tr(E£E?E?E*p) + £ [ Xmn Tr(E* *E?E? e£ P ) + X * mn Tr(E% E?E?E£p)]. 



m m<n 

Using E^E* = uf™E*E* and E^E* = w^E^E* this becomes: 



MT rsPk ) = Ti[ p kgip)] (Y. Lurm XmmTT(T rs p)+Y K"XmnTV(£ n At E*T rs p) + X * mn Tr(E^ E*T rs p)}) (7) 



rn<n 

r 



Therefore, we have an additional set of d equations to iden- that Tt(E^ E^T rsP ) ^ 0. Suppose we now measure another 
tify the off-diagonal elements of the superoperator, provided 
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unitary operator T r i s i = Ef,Ef, that commutes with Sy, i.e. 
T r i s i G N(S), and also commutes with T rs , and satisfies the 
relations T r , a , ^ T b s S lj a (where < a,b < d), Ef ^ I 
and Ef 7^ /. Such a measurement results in d equations 
for Tr(T r / s < ph), similar to those for Tr(T rs pfc). However, for 
these equations to be useful for characterization of the dynam- 
ics, one needs show that they are all linearly independent. In 
the next section, we find the maximum number of independent 
and commutating unitary operators T rs such that their expec- 
tation values on the output state, Tr(T rs pk), result in linearly 
independent equations to be d — 1, see Fig. [4] I.e., we find 
an optimal Abelian set of unitary operators such that the joint 
probability distribution functions of their eigenvalues and sta- 
bilizer eigenvalues at the output state are linearly independent. 



VI. LINEAR INDEPENDENCE AND OPTIMALITY OF 
MEASUREMENTS 

Before presenting the proof of linear independence of the 
functions Ti (T rs pf.) and of the optimality of the DCQD algo- 
rithm, we need to introduce the following lemmas and defini- 
tions. 

Lemma 3 If a stabilizer group, S, has a single generator, the 
order of its normalizer group, N(S), is d 3 . 

Proof: Let us consider the sole stabilizer generator S12 = 
EfEf, and a typical normalizer element T V2 > ~ EfE§, 
where Ef = X^Z?\ Ef = X^Z^, Ef, = X™ Z^ 1 and 
Ef = X q z' Z p z' . Since Sn and Ty 2 i commute, we have 

S 12 T V2 , = oj^P^'-^i'T V2 ,S 12 , where YLiPi^ " 
qip' v = (mod d). We note that for any particular code with a 
single stabilizer generator, all q\,p\,q 2 and p 2 are fixed. Now, 
by LemmaQ] for given values of q[,p'x and q' 2 there is only one 
value for p' 2 that satisfies the above equation. However, each 
of q[,Pi and q' 2 can have d different values. Therefore, there 
are d 3 different normalizer elements Tyy. 

Lemma 4 Each Abelian subgroup of a normalizer, which in- 
cludes the stabilizer group {Sfj} as a proper subgroup, has 
order d 2 . 

Proof: Suppose T rs is an element of N(S), i.e., it commutes 
with Sij . Moreover, all unitary operators of the form T b s Sij a , 
where < a,b < d, also commute. Therefore, any Abelian 
subgroup of the normalizer, A C N(S), which includes {Sfj } 
as a proper subgroup, is at least order of d 2 . Now let T r / S > 
be any other normalizer element, i.e., T r i s i ^ T b s Sfj with 
< a,b < d, which belongs to the same Abelian subgroup 
A. In this case, any operator of the form T b , s ,T b s Sfj would 
also belong to A. Then all elements of the normalizer should 
commute or A = N(S), which is unacceptable. Thus, either 
T r , s , = T^Sfj or T r , s , £ A, i.e., the order of the Abelian 
subgroup A is at most d 2 . 

Lemma 5 There are d + 1 Abelian subgroups, A, in the nor- 
malizer N (S). 



Proof: Suppose that the number of Abelian subgroups which 
includes the stabilizer group as a proper subgroup is n. Using 
Lemmas [3] and |H we have: d 3 = nd 2 — (n — i)d, where the 
term (n — l)d has been subtracted from the total number of 
elements of the normalizer due to the fact that the elements 
of the stabilizer group are common to all Abelian subgroups. 
Solving this equation for n, we find that n = d _{ = d+1. 

Lemma 6 The basis of eigenvectors defined by d + 1 Abelian 
subgroups of N(S) are mutually unbiased. 

Proof: It has been shown [30] that if a set of d 2 — 1 traceless 
and mutually orthogonal d x d unitary matrices can be par- 
titioned into d+1 subsets of equal size, such that the d — 1 
unitary operators in each subset commute, then the basis of 
eigenvectors corresponding to these subsets are mutually un- 
biased. We note that, based on Lemmas [3]S and [5] and in the 
code space (i.e., up to multiplication by the stabilizer elements 
{Sij}), the normalizer N(S) has d 2 ~ 1 nontrivial elements, 
and each Abelian subgroup A, has d — 1 nontrivial commut- 
ing operators. Thus, the bases of eigenvectors defined by d+ 1 
Abelian subgroups of N(S) are mutually unbiased. 

Lemma 7 Let C be a cyclic subgroup of A, i.e., C c A c 
N{S). Then, for any fixed T G A, the number of distinct left 
(right) cosets, TC (CT), in each A is d. 

Proof: We note that the order of any cyclic subgroup C C A, 

such as T^ s with < b < d, is d. Therefore, by Lemma|4] the 

,2 

number of distinct cosets in each A is = d. 

a 

Definition 2 We denote the cosets of an (invariant) cyclic sub- 
group, C a , of an Abelian subgroup of the normalizer, A v , by 
A v /C a , where v = 1,2, . . . ,d+l. We also represent generic 
members of A v /C a as T^ s Sfj , where < a,b < d. The mem- 
bers of a specific coset A v /C ao are denoted as T^Sff, where 
ao represents a fixed power of stabilizer generator Sij, that 
labels a particular coset A v /C aa , and b (0 < b < d) labels 
different members of that particular coset. 

Lemma 8 The elements of a coset, T^S"? (where T rs = 
EfEf, Sij = EfEf and < b < d) anticommute with 
Ef with different eigenvalues ui k . I.e., there are no two dif- 
ferent members of a coset, A v /C ao , that anticommute with 
Ef with the same eigenvalue. 

Proof: First we note that for each Tf. s = {Ef) b {Ef) b , 
the unitary operators acting only on the principal subsystem, 
(Ef) b , must satisfy either (a) (Ef) b = Ef or (b) (Ef) b ^ 
Ef. In the case (a), and due to [T rs ,Sij] = 0, we should 
also have (Ef ) b = Ef , which results in T b s = S tj ; i.e., 
T b s is a stabilizer and not a normalizer. This is unacceptable. 
In the case (b), in particular for b = 1, we have EfEf = 
cu ri EfEf. Therefore, for arbitrary b we have (Ef) b Ef = 
u br 'Ef(Ef ) b . Since < b < d, we conclude that u br ' ^ 
cu b Ti for any two different values of b and b' . 

As a consequence of this lemma, different (Ef) b , for < 
b < d, belong to different Vt^'s. 
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FIG. 4: A diagram of the error-detection measurement for estimating quantum dynamical coherence: we measure the sole stabilizer generator 
at the output state, by applying projection operators corresponding to its different eigenvalues P^. We also measure d — 1 commuting operators 
that belong to the normalizer group. Finally, we calculate the probability of each stabilizer outcome, and joint probability distributions of the 
normalizers and the stabilizer outcomes. Optimally, we can obtain d 2 linearly independent equations by appropriate selection of the normalizer 
operators as it is shown in the next section. 



Lemma 9 For any fixed unitary operator G W^ k , where 
k 7^ 0, and any other two independent operators E^ and E^ 
that belong to the same Wl , we always have uj Tm ^ oj Tn , 
where E^E^ = u rm E%E? and E*E£ = Lu r "E^Ef. 

Proof: We need to prove for operators E^,E^,E^ e W k 
(where k ^ 0), that we always have: E^ ^ E^ ==> ^ 
uj r ". Let us prove the converse: u) rm = Lo Tn => E^ = E^. 
We define Ef = X qi Z Pi , Ef = X q -Z Pr , E£ = X q ™Z Pm , 
E£ = X qn Z Pn . Based on the definition of subsets W l k with 
k ^ 0, we have: piq m - qip m = p t q n - q t p n = k (mod d) 
k ■ it! (I), where t is an integer number. We need to show if 

P r q m - qrPm = PrQn - QrPn = ^'(mod d) = k' + t' d (II), 

then Ei = E* 

We divide the equations (I) by qiq m or qiq n to get: |f = 
MM + Ehl = £kM + 2il (p). We also divide the equations (II) 
by^ m org r ^toget:"|- = ^ + ^ = ^ + ^ (IV). 
By subtracting the equation (IF) from (F) we get: q n ( k ~^ td — 
*L±^) = g m (^±M _ l^A) (1). Similarly, we can obtain 

the equation p n (^ - = Pm (^ - ^±f^) (2). 

Note that the expressions within the brackets in both equa- 
tions (1) or (2) cannot be simultaneously zero, because it will 
result in piq r — qip r = 0, which is unacceptable for k ^ 0. 



Therefore, the expression within the brackets in at least one of 
the equations (1) or (2) is non-zero. This results in q n = q m 
and/or p n = p m . Consequently, considering the equation (I), 
we have E^ x = E£. 



A. Linear independence of the joint distribution functions 

Theorem 1 The expectation values of normalizer elements 
on a post-measurement state, pk, are linearly independent if 
these elements are the d — 1 nontrivial members of a coset 
A v /C ao . I.e., for two independent operators T rs , T r > s > G 
A v /C aa , we have Tr^TrsPk) ^ c Tr(T r i s > p^), where c is an 
arbitrary complex number. 

Proof: We know that the elements of a coset can be written 
as T b s S% = {EfEffSfi, where b = 1, 2, . . . , d - 1. We 
also proved that {Ef) b belongs to different W\. (k ^ 0) for 
different values of b (see Lemma[8]). Therefore, according to 
Lemma [9] and regardless of the outcome of k (after measur- 
ing the stabilizer 5V, ), there exists one member in the coset 
A v /C ao that has different eigenvalues u> rm with all (indepen- 
dent) members G Wj,. The expectation value of T b s S^° 
is: 



Tv(T b s S«° Pk ) = ]T Xmm T,{EfjT b rs St°Eip), + £ [ Xmn Tr(E^T* s S?? E £p) + X * mn Tr(E^T b s S% E £ P )\, (8) 

m m<n 

Tr(T r b sPfc ) = J2 u^Xmm MT b rsP ) + £ l^Xmn Tr(E? EiT b rs p) + J^ X * mn Tr(E^ E*T b sP )}, (9) 



where Lj Vm ^ uj t ' 1 ^= ... for all elements E^,E^, . . . that 
belong to a specific W\. Therefore, for two independent 
members of a coset denoted by b and b' (i.e., b ^ b'), we 
have (uj b ' rm , Lu b ' r " , . . .) ^ c [uo brm , w br " , . . .) for all values 
of < b,b' < d, and any complex number c. We also note 
that we have Tr(E^E^ n T b s p) ? c Tr E*T b ' sP ), since 
T b s ~ b is a normalizer, not a stabilizer element, and its action 
on the state cannot be expressed as a global phase. Thus, for 



any two independent members of a coset A v /C ao , we always 
haveTr^p^cTrtT^). 

In summary, after the action of the unknown dynamical 
process, we measure the eigenvalues of the stabilizer gener- 
ator, EfE^, that has d eigenvalues for k = 0, 1, . . . , d — 1 
and provides d linearly independent equations for the real and 
imaginary parts of Xmn- This is due to that the outcomes 
corresponding to different eigenvalues of a unitary operator 
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are independent. We also measure expectation values of all 
the d — 1 independent and commuting normalizer operators 
Tr S Sij € A v /C ao , on the post-measurement state p k , which 
provides (d — 1) linearly independent equations for each out- 
come k of the stabilizer measurements. Overall, we obtain 
d + d(d — 1) = d 2 linearly independent equations for charac- 
terization of the real and imaginary parts of Xmn by a single 
ensemble measurement. In the following, we show that the 
above algorithm is optimal. I.e., there does not exist any other 
possible strategy that can provide more than log 2 d 2 bits of in- 
formation by a single measurement on the output state £(p). 

B. Optimality 

Theorem 2 The maximum number of commuting normalizer 
elements that can be measured simultaneously to provide lin- 



ear independent equations for the joint distribution functions 
Tr(T b s S«p k )isd-l. 



Proof: Any Abelian subgroup of the normalizer has order 
d 2 (see Lemma 0). Therefore, the desired normalizer opera- 
tors should all belong to a particular A v and are limited to d 2 
members. We already showed that the outcomes of measure- 
ments for d — 1 elements of a coset A v /C a , represented by 
T^Sfj (with b ^ 0), are independent (see TheoremQ]!. Now 

we show that measuring any other operator, T b s S"j , from any 
other coset A v /C a ', results in linearly dependent equations 
for the functions w =tr(T r b s 5'° / 9fe) and w' =tr(T b s S^ p k ) as 
the following: 



w = Tr(T? s S: jPk ) = £ Xmm ME^T r b s S*E*p) + £ [ Xmn ^{E^T^E^p) + X * mn Tr(E^T b s S^ p)] 

m m<n 

w' = Tr(T? s S?;p k ) = E Xm m Tr(E^T b s S^Eip) + £ [ Xmn Tr{E^T^E^p) + X * mn tv{E^T b rs S?' 3 E*p)]. 

m m<n 

Using the commutation relations T b s S^E^ = u br ™ +ai ™ E^T b s S^,we obtain: 

w = J2" brm+mm Xmm Tr(T* s p) + E [u br ~ +m ~Xmn Tr(E^ E? n T b s p) + w^+ ai « X * mn Tr(E^ E A T? s p)} 
w' = Y J U bTm+a ' lm Xmm tr{T b rs p) + E [u brm+a ' tm Xmn Tt{E$ t E f n T b s p) + x * mn Tr(E^ Ei T b rs p)\ , 

m m<n 

where we also used the fact that both Sfj and Sf^ are stabilizer elements. Since all of the operators E A belong to the same W k , 
we have i m = i n = k, and obtain: 

W = u; ak ( J2^ rm Xmm Tv(T b s p) + J2 \^ rm Xmn Tv(E^ E^sP) + U,*" X * mn Tr(E$ E£T? sP )] J 

w ' = u a ' k y£uj br - Xmm Tv{T b sP ) + £ [^Xmn Tr{E^EiT b s p)+uj br - X * mn Tr(E^E^ s p)) J . 

i 



Thus, we have w' = uj^ a ~ a *> k w, and consequently the 
measurements of operators from other cosets A v jC a > do 
not provide any new information about Xmn beyond the 
corresponding measurements from the coset A v /C a . 

For another proof of the optimality, based on fundamen- 
tal limitation of transferring information between two parties 
given by the Holevo bound see Ref. $2M- In principle, one can 
construct a set of non- Abelian normalizer measurements, from 
different A v , where w = l,2,...,d+l,to obtain information 
about the off-diagonal elements Xmn . However, determining 
the eigenvalues of a set of noncommuting operators cannot 
be done via a single measurement. Moreover, as mentioned 



above, by measuring the stabilizer and d — 1 Abelian nor- 
malizes, one can obtain log 2 d 2 bits of classical information, 
which is the maximum allowed by the Holevo bound l3lll . 
Therefore, other strategies involving non-Abelian, or a mix- 
ture of Abelian and non-Abelian normalizer measurements, 
cannot improve the above scheme. It should be noted that 
there are several possible alternative sets of Abelian normal- 
izers that are equivalent for this task, we address this issue in 
the next lemma. 

Lemma 10 The number of alternative sets of Abelian normal- 
izer measurements that can provide optimal information about 
quantum dynamics, in one ensemble measurement, is d 2 . 
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Proof: We have d + 1 Abelian normalizers A v (see Lemma 
[5]). However, there are d of them that contain unitary opera- 
tors that act nontrivially on both qudit systems A and B, i.e., 
T 7 b s = {E A Eff, where E A ^ I, Ef ^ 7. Moreover, in 
each we have <i cosets (see Lemma|5]) that can be used for 
optimal characterization of Xmn- Overall, we have d 2 pos- 
sible sets of Abelian normalizers that are equivalent for our 
purpose. 

In the next section, we develop the algorithm further to ob- 
tain complete information about the off-diagonal elements of 
the superoperator by repeating the above scheme for different 
input states. 



VII. REPEATING THE ALGORITHM FOR OTHER 
STABILIZER STATES 

we have shown that by performing one ensemble measure- 
ment one can obtain d 2 linearly independent equations for 
Xmn- However, a complete characterization of quantum dy- 
namics requires obtaining d 4 — d 2 independent real parameters 
of the superoperator (or d 4 for non-trace preserving maps), 
we next show how one can obtain complete information by 
appropriately rotating the input state and repeating the above 
algorithm for a complete set of rotations. 

Lemma 11 The number of independent eigenkets for the er- 
ror operator basis {Ej}, where j — 1, 2, . . . , d 2 — 1, is d+ 1. 
These eigenkets are mutually unbiased. 

Proof: We have d 2 — 1 unitary operators ,Ei. We note that 
the operators Ef for all values of 1 < a < e? — 1 commute and 
have a common eigenket. Therefore, overall we have (d 2 — 
l)/(d — 1) = d + 1 independent eigenkets. Moreover, it has 
been shown lf30h that if a set of d 2 — 1 traceless and mutually 
orthogonal dxd unitary matrices can be partitioned into d+ 1 
subsets of equal size, such that the d — 1 unitary operators in 
each subset commute, then the basis of eigenvectors defined 
by these subsets are mutually unbiased. 

Let us construct a set of d + 1 stabilizer operators E A E B , 
such that the following conditions hold: (a) E A ,E B ^ I, 

(b) (E A ) a ^ E A for i ^ i' and 1 < a < d - 1. Then, 
by preparing the eigenstates of these d + 1 independent sta- 
bilizer operators, one at a time, and measuring the eigen- 
values of Sij and its corresponding d — 1 normalizer oper- 
ators T^ s Sfj £ A v /C a , one can obtain (d + l)d 2 linearly 
independent equations to characterize the superoperator's off- 
diagonal elements. The linear independence of these equa- 
tions can be understood by noting that the eigenstates of all 
operators Ef of these d + 1 stabilizer operator Sij are mutu- 
ally unbiased (i.e., the measurements in these mutual unbiased 
bases are maximally non-commuting). For example the bases 

{|0>,|1>}, {{\+)x>\-)x} and {|+)y.|-M eigenstates 
of the Pauli operators Z, X, and Y) are mutually unbiased, 
i.e., the inner products of each pair of elements in these bases 
have the same magnitude. Then measurements in these bases 
are maximally non-commuting ll32ll . 



To obtain complete information about the quantum dynam- 
ical coherence, we again prepare the eigenkets of the above 
d + 1 stabilizer operators E A E B , but after the stabilizer mea- 
surement we calculate the expectation values of the operators 
T^, 8 ,Sfj belonging to other Abelian subgroups A v i )C a of the 
normalizer, i.e., A v > ^ A v . According to Lemma|6]the bases 
of different Abelian subgroups of the normalizer are mutu- 
ally unbiased, therefore, the expectation values of T^, s , Sf^ and 
T^Sfj from different Abelian subgroups A v > and A v are in- 
dependent. In order to make the stabilizer measurements also 
independent we choose a different superposition of logical ba- 
sis in the preparation of d + 1 possible stabilizer state in each 
run. Therefore in each of these measurements we can obtain 
at most d 2 linearly independent equations. By repeating these 
measurements for d — 1 different A v over all d + 1 possible 
input stabilizer state, we obtain (d + l)(d — l)d 2 = d 4 — d 2 
linearly independent equations, which suffice to fully charac- 
terize all independent parameters of the superoperator's off- 
diagonal elements. In the next section, we address the general 
properties of these d + 1 stabilizer states. 



VIII. GENERAL CONSTRAINTS ON THE STABILIZER 
STATES 

The restrictions on the stabilizer states p can be expressed 
as follows: 

Condition 1 The state p = \4>ij) (4>ij | is a non-separable pure 
state in the Hilbert space of the two-qudit system H. I.e., 

Condition 2 The state | <fiij ) is a stabilizer state with a sole sta- 
bilizer generator Sij = EjEj. I.e., it satisfies S*°- \4>ij) = 
u> ak \4>ij), where k G {0, 1, . . . , d — 1} denotes a fixed eigen- 
value of Sij, and a = 1, . . . , d — 1 labels d — 1 nontrivial 
members of the stabilizer group. 

The second condition specifies the stabilizer subspace, Vs, 
that the state p lives in, which is the subspace fixed by all 
the elements of the stabilizer group with a fixed eigenvalues 
k. More specifically, an arbitrary state in the entire Hilbert 

d-l 

space H can be written as \<fi) = a uu' \ u )a W) b 

u.u' — 

where {\u)} and {|it')} are bases for the Hilbert spaces 
of the qudits A and B, such that X q \u) = \u + q) and 
Z p | it) = uj pu \u). However, we can expand \<fr) in another 

d-l 

basis as |</>) = fiw 1 \ v ) a \ v ')b' sucn m at X q \v) = 

v^v'—Q 

oj qv \v) and Z p \v) = \v +p). Let us consider a stabilizer 
state fixed under the action of a unitary operator EfE? = 

{X A Y{X B ) q (Z a )p(Z b )p' with eigenvalue co k . Regardless 
of the basis chosen to expand \(f>ij), we should always have 
Sij \4>ij) = to k \4>ij)- Consequently, we have the constraints 
pu (B p'u' = k, for the stabilizer subspace Vs spanned by the 
{\u) <g> \u'}} basis, and q(v © p) © q'(v' © p') = k, if V s is 
spanned by {\v) © basis, where © is addition mod(<i). 
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From these relations, and also using the fact that the bases 
{\v)} and {\u}} are related by a unitary transformation, one 
can find the general properties of Vs for a given stabilizer gen- 
erator E^Ej and a given k. 

We have already shown that the stabilizer states p should 
also satisfy the set of conditions Tx[E^ E^p] ^ and 
Tv{E^E^T b s p) ^ for all operators belonging to the 
same W%., where T b s (0 < b < d — 1) are the members 
of a particular coset A v /C a of an Abelian subgroup, A v , of 
the normalizer N(S). These relations can be expressed more 
compactly as: 

Condition 3 For stabilizer state p — \4>ij) (<fiij\ = \4>c) (4>c\ 
and for all G W l k we have: 

(<f> c \E^E£T b s \J> c )^0, (10) 
where here < b < d — 1 . 

Before developing the implications of the above formula 
for the stabilizer states we give the following definition and 
lemma. 

Definition 3 Let be the logical basis of the code space 

that is fixed by the stabilizer generator E^E^. The stabilizer 

d-1 

state in that basis can be written as \(j> c ) — a i I0i> an d a ^ 

i=o 

the normalizer operators, T rs , can be generated from tensor 
products of logical operations X and Z defined as Z \ l) L = 
J \l) L <mdX\l) L = + For example: \l) L = \k)\k), 
~Z = Z®Iax\&X~=X®X, where X \k) = \k + 1) and 

Z\k) =LO k \k). 

Lemma 12 For a stabilizer generator E^E^ and all unitary 
operators E^ G W k , we always have E^E^ = ui c Z a , 
where Z is the logical Z operation acting on the code space 
and a and c are integers. 

Proof: Let us consider Ef = X qi Z Pi , and two generic op- 
erators E£ and E£ that belong to W l k : E% = X gm Z Pm and 
E£ = Xi" ZP" . From the definition of W l k (see Definition[D 
we have p t q m - q t p m = piq n - qip n = k (modd) = k + td. 
We can solve these two equations to get: q m —<ln = qiiPmQn— 
q m Pn)/(k + td) andp m -p n = Pi(p m q n - q m p n ) / (k + td) . 
We also define p m q n — q m p n = k' + t'd. Therefore, we obtain 
Qm — Qn = 1i a an d Pm —Pn= Pi a , where we have introduced 

a= (k' + t'd)/(k + td). (11) 

Moreover, we have E^ = X^" d ~ q ^ Z^" d ^ p "^ for some 
other integer t". Then we get: 

E^E* = uj c X( t " d +i™-<* s >Z( t " d+ P™-^ 

= ^x^-^z^-^ 
= u} c (x q *z p *y\ 

where c= (t"d-p n )(t"d+q m - q n ). However, X^Z P ' ®I 



acts as logical Z on the code subspace, which is the eigenstate 

of EfEf. Thus, we obtain E^ E^ = uj c Z a . 

Based on the above lemma, for the case of b = we obtain 

d-1 

(0c| E^Ei \cj> c ) = cj c (<f> c \Z a \d> c ) =^Y. " al \ai | 2 . 

;=o 

Therefore, our constraint in this case becomes 
Sfc=o u>al \ a i | 2 7^ 0' which is not satisfied if the stabi- 
lizer state is maximally entangled. For b ^ 0, we note that 
T b s are in fact the normalizers. By considering the general 
form of the normalizer elements as T b s = (X Z ) b , where q, 
p G {0, 1, . . . , d — 1}, we obtain: 

{M E^EiT b rs \d> c ) = c c {<j> c \ T(X q Z P ) b fa) 

k=0 

d-1 
1=0 

Overall, the constraints on the stabilizer state, due to condition 
(iii), can be summarized as: 

d-1 

Y,^ a+bp)l a:a l+bq ^0 (12) 

1=0 

This inequality should hold for all b G {0,1, . . . ,d — 1}, 
and all a defined by Eq. (fTTT i. however, for a particular coset 
A v /C a the values of q and p are fixed. One important prop- 
erty of the stabilizer code, implied by the above formula with 
b = 0, is that it should always be a nonmaximally entangled 
state. In the next section, by utilizing the quantum Hamming 
bound, we show that the minimum number of physical qudits, 
n, needed for encoding the required stabilizer state is in fact 
two. 



IX. MINIMUM NUMBER OF REQUIRED PHYSICAL 
QUDITS 

In order to characterize off-diagonal elements of a superop- 
erator we have to use degenerate stabilizer codes, in order to 
preserve the coherence between operator basis elements. De- 
generate stabilizer codes do not have a classical analog JH]. 
Due to this fact, the classical techniques used to prove bounds 
for non-degenerate error-correcting codes cannot be applied 
to degenerate codes. In general, it is yet unknown if there are 
degenerate codes that exceed the quantum Hamming bound 
|Q]]. However, due to the simplicity of the stabilizer codes 
used in the DCQD algorithm and their symmetry, it is possi- 
ble to generalize the quantum Hamming bound for them. Let 
us consider a stabilizer code that is used for encoding k logi- 
cal qudits into n physical qudits such that we can correct any 
subset of t or fewer errors on any n e ^ n of the physical qu- 
dits. Suppose that ^ j ^ t errors occur. Therefore, there are 
( n f ) possible locations, and in each location there are (d 2 — 1 ) 
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different operator basis elements that can act as errors. The to- 
tal possible number of errors is E*=o (™T ) — l) 5 - ^ me 
stabilizer code is non-degenerate each of these errors should 
correspond to an orthogonal d k -dimensional subspace; but if 
the code is uniformly g-fold degenerate (i.e., with respect to 
all possible errors), then each set of g errors can be fit into 
an orthogonal d k -dimensional subspace. All these subspaces 
must be fit into the entire d" -dimensional Hilbert space. This 
leads to the following inequality: 



t 

E 

3=0 



(d 2 - l)id k 



(13) 



We are always interested in finding the errors on one phys- 
ical qudit. Therefore, we have n e = 1, j 6 {0,1} and 

( c " c ) = 1, and Eq. <T3J becomes Y? j=0 < d n . 

In order to characterize diagonal elements, we use a nonde- 
generate stabilizer code with n = 2, k = and g = 1, 
and we have ^2j =0 (d 2 — l) 3 = d 2 . For off-diagonal ele- 
ments, we use a degenerate stabilizer code with n = 2, k = 1 

and g = d, and we have Ej=o ~ ~tr ~ = Therefore, 
in the both cases the upper-bound of the quantum Hamming 
bound is satisfied by our codes. Note that if instead we use 
n = k, i.e., if we encode n logical qudits into n separa- 
ble physical qubits, we get Ej=o g — T This can 
only be satisfied if g = d 2 , in which case we cannot ob- 
tain any information about the errors. The above argument 
justifies Condition (i) of the stabilizer state being nonsepara- 
ble. Specifically, it explains why alternative encodings such 
as n = k = 2 and n = k = 1 are excluded from our 
discussions. However, if we encode zero logical qubits into 
one physical qubit, i.e., n = 1, k = 0, then, by using a 

d-fold degenerate code, we can obtain E,J=o ~ = ^ 

which satisfies the quantum Hamming bound and could be 
useful for characterizing off-diagonal elements. For this to 
be true, the code \<f> c ) should also satisfy the set of condi- 
tions (<f> c \E^E A \<j> c ) ? and (0 C | E A * E*T$ a \<j> D ) f 0. 
Due to the d-fold degeneracy of the code, the condition 
((j> c \ E A ^ E^ \<j) c ) ^ is automatically satisfied. However, 
the condition (0 C | E A ^E A T h rs \<j> c ) ^ can never be satisfied, 
since the code space is one-dimensional, i.e., d k = 1, and the 
normalizer operators cannot be defined. I.e., there does not 
exist any nontrivial unitary operator T^ s that can perform log- 



ical operations on the one-dimensional code space. 

we have demonstrated how we can characterize quantum 
dynamics using the most general form of the relevant stabi- 
lizer states and generators. In the next section, we choose a 
standard form of stabilizers, in order to simplify the algorithm 
and to derive a standard form of the normalizer. 



X. STANDARD FORM OF STABILIZER AND 
NORMALIZER OPERATORS 

Let us choose the set {|0) , |1) , \k — 1)} as a standard 



basis, such that Z\k) 



\k) and X \k) 



1 



In 



order to characterize the quantum dynamical population, we 
choose the standard stabilizer generators to be {X A X B ) q 
and {Z A (Z B ) d ~ 1 ] p . Therefore, the maximally entangled in- 

d-l 

put states can be written as \ip c ) = -4= ^ jfc)^ \k) B . In 

v k=0 

order to characterize the quantum dynamical coherence we 
choose the sole stabilizer operator as [Ef(E B ) d ~ 1 ] a , which 

d-l 



E a i V)a where 

i=0 

Ei \i) = oj 1 \i) and \i) represents one of d + 1 mutually unbi- 
ased basis states in the Hilbert space of one qudit. The normal- 
izer elements can be written as T b qp = {X q Z P ) b E A v JC ao , 
for all < b < d - 1, where X = Ei® Ea ,~Z = Ei® I, 



has an eigenket of the form \tp c ) 



Ei 



\i + 1) and Ei \i) 



and A Vo /C ao represents 



a fixed coset of a particular Abelian subgroup, A VQ , of the nor- 
malizer N(S). For example, for a stabilizer generator of the 
form [Ef(E?) d - 1 ] a = [Z A (Z B ) d - 1 }P we prepare its eigen- 

d-l 

ket \ f c ) = E a k W) a W)b i anc ^ tne normalizers become 

k=Q 

T h qv = (X q Z P ) b , where X = X ® X and Z = Z ® I. Using 
this notations for stabilizer and the normalizer operators, we 
provide an overall outline for the DCQD algorithm in the next 
section. 



XI. ALGORITHM: DIRECT CHARACTERIZATION OF 
QUANTUM DYNAMICS 

The DCQD algorithm for the case of a qudit system is sum- 
marized as follows (see also Figs.|5]and|6]): 



J 



,, n E A pE A \ where p denotes the state of the 



Inputs: (1) An ensemble of two-qudit systems, A and B, prepared in the state \0) A ® \0) B . (2) An arbitrary unknown CP 

quantum dynamical map E, whose action can be expressed by E(p) = Em n=o X- 
primary system and the ancilla. 

Output: £, given by a set of measurement outcomes in the procedures (a) and (b) below: 
Procedure(a): Characterization of quantum dynamical population (diagonal elements Xmm of x)< see Fig. [5] 
1. Prepare \ipo) = \0) A ® \0) B , a pure initial state. 

d-l 



2. Transform it to \ip c ) 



v^7 



E \k) A \k) B , a maximally entangled state of the two qudits. This state has the stabilizer 



fc=0 



operators E A E B = (X A X B ) q and E A E B , = [Z A (Z B ) d - 1 ] p for 0<p,q<d-l. 



3. Apply the unknown quantum dynamical map to the qudit A: £(p) 



2^,m,n=0 Xri 



E A pE A \ where p = \4> c ) (4> c 
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FIG. 5: Procedure (a): Measuring the quantum dynamical population (diagonal elements Xmm). The arrows indicate direction of time. (1) 

d-l 

Prepare |^o) = |0)a ® |0)b> a P ure initial state. (2) Transform it to \tpc) — -75 ~}2 \k) A \k) g , a maximally entangled state of the two qudits. 

V fc=0 

This state has the stabilizer operators S = X A X B and S' — Z (Z B ) . (3) Apply the unknown quantum dynamical map to the qudit 
A, £(p), where p — \4>c) {4>c\- (4) Perform a projective measurement PkPk' , for k, k' = 0, . . . , d — 1, corresponding to eigenvalues of 
the stabilizer operators S and S' . Then calculate the joint probability distributions of the outcomes k and k': Tr[PkPk'£(p)] = Xmm- The 
elements Xmm represent the population of error operators that anticommute with the stabilizer generators S and S' with eigenvalues uj k and 
uj k , respectively. The number of ensemble measurements for procedure (a) is one. 

4. Perform a projective measurement P k P k , : £{p) h-> P k P k ,£(p)P k P k ,, where P k = ± ^f=o u~ lk (EfEf) 1 , and P k , = 
2 2c=o 10 ~ l k (EftE^) 1 , and calculate the joint probability distributions of the outcomes k and k': 

Tr{P k P k ,£(p)} 

Xmm ■ 

Number of ensemble measurements for Procedure (a): 1. 

Procedure (b): Characterization of quantum dynamical coherence (off-diagonal elements Xmn of x)< see Fig- 

1. Prepare \<po) = \0) A ® |0) B , a pure initial state. 

d-l 

2. Transform it to \cp c ) = ^ on \i) A \i) B , a non-maximally entangled state of the two qudits. This state has stabilizer 

i=0 

operators [Ef(E? ) d " 1 ] a . 

3. Apply the unknown quantum dynamical map to the qudit A: £{p) = ^~^ Xmn E^pE^ , where p = \(f> c ) ((f> c \. 

4. Perform a projective measurement 

P k : £{p) -> Pk = P k £( P )P k = Xmm E&pE% + ]T ( Xmn ^pS^t + x * mn E*pE%), 

m rn < n 

where P fc = 3 E^tii u~ lk (EfEf) 1 and £^ = X* m Z Pm £ W{, and calculate the probability of outcome k: 

Tr[P k £(p)} = Xmm + 2 MXmn Tr(E^ E^p)}. (14) 

m m<n 

5. Measure the expectation values of the normalizer operators T^ p = (X q Z P ) b £ A Va /C aa , for all < 6 < <i — 1, where 

X = Ei (g> Ei , Z = Ei <E> I , Ei \i) = u> % \i), Ei \i) = \i+ 1), where A V(l /C ao represents a fixed coset of a particular Abelian 
subgroup, A vo , of the normalizer N(S). 

Tr(T g » = Y^^Xmm Tr(T r b sP ) + £ [co^^Xmn Hr(E? Effi.p) + ^- w "xL Tr{E^ E*T b rs p)}. 

m m<n 

6. Repeat the steps (l)-(5) d + 1 times, by preparing the eigenkets of other stabilizer operator [Ef(E^) d ^ 1 ] a for all i S 
{1, 2, d + 1}, such that states |i) B in the step (2) belong to a mutually unbiased basis. 

7. Repeat the step (6) up to d — 1 times, each time choosing normalizer elements T b p from a different Abelian subgroup 
A v /C a , such that these measurements become maximally non-commuting. 

Number of ensemble measurements for Procedure (b): (d + l)(d — 1). 
Overall number of ensemble measurements: d 2 . 

Note that at the end of each measurement in Figs. Bland Ip, 1 , ., , . . , . £ .... 

r-r, , . , , i v ) \i) a W r - is the common eigenket of the stabilizer gen- 

\b\ the output state - a maximally entangled state, \<pe) = £r 
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FIG. 6: Procedure (b): Measuring the quantum dynamical coherence (off-diagonal elements Xmn)- (1) Prepare | t^o ) = |0} A ® |0) fl , a pure 

d-l 

initial state. (2) Transform it to \ipc) = ai \i) a l*)s> a nonmaximally entangled state of the two qudits. This state has a sole stabilizer 

j=0 

operator of the form S = Ef'{E^ 3 ) d ~ 1 . (3) Apply the unknown quantum dynamical map to the qudit A, £(p), where p = \<j>c) (<j>c\- 
(4) Perform a projective measurement P^, and calculate the probability of the outcome k: tr[Pk£(p)]. (5) Measure the expectation values of 
all normalizer operators T rs that simultaneously commute with the stabilizer generator S. There are only d — 1 such operators T rs that are 
independent of each other, within a multiplication by a stabilizer generator; and they belong to an Abelian subgroup of the normalizer group. 
(6) Repeat the steps (l)-(5) d + 1 times, by preparing the eigenkets of other stabilizer operator Ef(Ef) d ^ 1 for alH € {1, 2, . . . , d + 1}, such 
that states \i) A in step (2) belong to a mutually unbiased basis l32ll . (7) Repeat the step (6) up to d — 1 times, each time choosing normalizer 
elements T rs from a different Abelian subgroup of the normalizer, such that these measurements become maximally non-commuting, i.e., their 
eigenstates form a set of mutually unbiased bases. The number of ensemble measurements for Procedure (b) is (d + l)(d — 1). 



erator and its commuting normalizer operators. For the pro- 
cedure (a), this state can be directly used for other mea- 
surements. This is indicated by the dashed lines in Fig. [5] 
For the procedure (b), the state \<Pe) can be unitarily trans- 
formed to another member of the same input stabilizer code, 

d-l 

\<pc) = a i I i) a I *)s' before another measurement. There- 
in 

fore, all the required ensemble measurements, for measuring 
the expectation values of the stabilizer and normalizer opera- 
tors, can always be performed in a temporal sequence on the 
same pair of qudits. 

In the previous sections, we have explicitly shown how the 
DCQD algorithm can be developed for qudit systems when d 
is prime. In the appendix [A] we demonstrate that the DCQD 
algorithm can be generalized to other A^-dimensional quan- 
tum systems with N being a power of a prime. 



XII. SUMMARY 

For convenience, we provide a summary of the DCQD al- 
gorithm. The DCQD algorithm for a qudit, with d being a 
prime, was developed by utilizing the concept of an error 
operator basis. An arbitrary operator acting on a qudit can 
be expanded over an orthonormal and unitary operator basis 



dS, 



{E^Ex,. . .,E d 2_ 1 }, where E = I and tr{EjEj) 
Any element £?, can be generated from tensor products of X 
and Z, where X \k) = \k + 1) and Z \k) = uj k \k), such that 
the relation XZ = 



ZX is satisfied 02811 . Here w is a dth 
root of unity and X and Z are the generalizations of Pauli 
operators to higher dimension. 

Characterization of Dynamical Population - A measure- 
ment scheme for determining the quantum dynamical popu- 
lation, Xmm> m a single experimental configuration. Let us 
prepare a maximally entangled state of the two qudits \<pc) — 



-4= \k)^\k) B . This state is stabilized under the action of 

stabilizer operators S = X A X B and S'' = Z A (Z B ) d ~ 1 , and 
it is referred to as a stabilizer state 01 13. After applying 
the quantum map to the qudit A, £(p), where p = \<f>c) (<^c|, 
we can perform a projective measurement PkPk'£(p)PkPk', 



d-l , ,-lk cl 
1=0 



S l ,P k ' = ^J2trJ co- l ' k 'S n ',mid 



where Pk 

lu = e aix l d . Then, we calculate the joint probability distri- 
butions of the outcomes k and k'\ Tr[PkPk'£(p)] = Xmm, 
where the elements Xmm represent the population of error op- 
erators that anticommute with stabilizer generators S and S' 
with eigenvalues u k and u k , respectively. Therefore, with a 
single experimental configuration we can identify all diagonal 
elements of superoperator. 

Characterization of Dynamical Coherence- For measuring 
the quantum dynamical coherence, we create a nonmaximally 

d-l 

entangled state of the two qudits \<pc) = S a i V) a V)b- 

i=0 

This state has the sole stabilizer operator S = Ef (E^) d ~ 1 
(for detailed restrictions on the coefficients on see Sec. I Villi ). 
After applying the dynamical map to the qudit A, £(p), we 
perform a projective measurement pu = Pk£(p)Pk, and 
calculate the probability of the outcome k: Tr[Pk£(p)] = 
J2 m X mm + 2 £ m< „ Rc[ Xmn Tr(E^E*p)}; where are 
all the operators in the operator basis, {E A }, that anticom- 
mute with the operator Ef with the same eigenvalue u k . We 
also measure the expectation values of all independent op- 
erators T rs = E A Ef of the Pauli group (where E A ^ I; 
Ef 7^ /) that simultaneously commute with the stabilizer 
generator S: Ti(T rs pk). There are only d — 1 such operators 
T rs that are independent of each other, within a multiplication 
by a stabilizer generator; and they belong to an Abelian sub- 
group of the normalizer group. The normalizer group is the 
group of unitary operators that preserve the stabilizer group by 



15 



conjugation, i.e., TST^ = S. We repeat this procedure d + 1 
times, by preparing the eigenkets of other stabilizer operator 
Ef-(Ef) d ~ 1 for all i G {1, 2, . . . , d+1}, such that states \i) A 
in input states belong to a mutually unbiased basis 13211 . Also, 
we can change the measurement basis d — 1 times, each time 
choosing normalizer elements T rs from a different Abelian 
subgroup of the normalizer, such that their eigenstates form 
a mutually unbiased basis in the code space. Therefore, we 
can completely characterize quantum dynamical coherence by 
(d+ l)(rf— 1) different measurements, and the overall number 
number of experimental configuration for a qudit becomes d 2 . 
For A^-dimensional quantum systems, with N a power of a 
prime, the required measurements are simply the tensor prod- 
uct of the corresponding measurements on individual qudits - 
see AppendixlAl For quantum system whose dimension is not 
a power of a prime, the task can be accomplished by embed- 
ding the system in a larger Hilbert space whose dimension is 
a prime. 

XIII. OUTLOOK 

An important and promising advantage of DCQD is for use 
in partial characterization of quantum dynamics, where in 
general, one cannot afford or does not need to carry out a full 
characterization of the quantum system under study, or when 
one has some a priori knowledge about the dynamics. Using 
indirect methods of QPT in those situations is inefficient, be- 
cause one has to apply the whole machinery of the scheme to 
obtain the relevant information about the system. On the other 
hand, the DCQD scheme has built-in applicability to the task 
of partial characterization of quantum dynamics. In general, 
one can substantially reduce the overall number of measure- 
ments, when estimating the coherence elements of the super- 
operator for only specific subsets of the operator basis and/or 
subsystems of interest. This fact has been demonstrated in 
Ref. lf26ll in a generic fashion, and several examples of partial 
characterization have also been presented. Specifically, it was 
shown that DCQD can be efficiently applied to (single- and 
two-qubit) Hamiltonian identification tasks. Moreover, it is 
demonstrated that the DCQD algorithm enables the simulta- 



neous determination of coarse-grained (semiclassical) physi- 
cal quantities, such as the longitudinal relaxation time Ti and 
the transversal relaxation (or dephasing) time T2 for a single 
qubit undergoing a general CP quantum map. The DCQD 
scheme can also be used for performing generalized quan- 
tum dense coding tasks. Other implications and applications 
of DCQD for partial QPT remain to be investigated and ex- 
plored. 

An alternative representation of the DCQD scheme for 
higher-dimensional quantum systems, based on generalized 
Bell-state measurements will be presented in Ref. 13311 - The 
connection of Bell-state measurements to stabilizer and nor- 
malizer measurements in DCQD for two-level systems, can be 
easily observed from Table II of Ref. yfl. Our presentation of 
the DCQD algorithm assumes ideal (i.e., error-free) quantum 
state preparation, measurement, and ancilla channels. How- 
ever, these assumptions can all be relaxed in certain situations, 
in particular when the imperfections are already known. A 
discussion of these issues is beyond the scope of this work 
and will be the subject of a future publication lf33ll . 

There are a number of other directions in which the results 
presented here can be extended. One can combine the DCQD 
algorithm with the method of maximum likelihood estimation 
113511 , in order to minimize the statistical errors in each experi- 
mental configuration invoked in this scheme. Moreover, a new 
scheme for continuous characterization of quantum dynamics 
can be introduced, by utilizing weak measurements for the re- 
quired quantum error detections in DCQD ll3ol l37ll . Finally, 
the general techniques developed for direct characterization 
of quantum dynamics could be further utilized for control of 
open quantum systems lf38tl . 
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APPENDIX A: GENERALIZATION TO ARBITRARY OPEN QUANTUM SYSTEMS 



Here, we first demonstrate that the overall measurements for a full characterization of the dynamics of an n qudit systems 
(with d being a prime) become the tensor product of the required measurements on individual qudits. One of the important 
examples of such systems is a QIP unit with r qubits, thus having a 2 '"-dimensional Hilbert space. Let us consider a quantum 
system consisting of r qudits, p = pi g) p2 <8> • ■ ■ <& p r , with a Hilbert space of dimension N ~ d r . The output state of such a 

system after a dynamical map becomes £(p) = Y^,m n=o Xmn EmpE^ where here {E m } are the unitary operator basis elements 
of an iV-dimensional Hilbert space. These unitary operator basis elements can be written as E m = X qm ^ Z pm i <E> X qm ^ Z Pm ? ® 
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TABLE I: Required physical resources for the QPT schemes: Standard Quantum Process Tomography (SQPT), Ancilla-Assisted Process To- 
mography using separable joint measurements (AAPT), using mutual unbiased bases measurements (MUB), using generalized measurements 
(POVM), see Ref. yfl, and Direct Characterization of Quantum Dynamics (DCQD). The overall number of measurements is reduced quadrat- 
ically in the DCQD algorithm with respect to the separable methods of QPT. This comes at the expense of requiring entangled input states, 
and two-qudit measurements of the output states. The non-separable AAPT schemes require many-body interactions that are not available 
experimentally (3]. 



Scheme 



dim(-H)° 



A', 



measurements 



required interactions 



SQPT 
AAPT 

AAPT (MUB) 
AAPT (POVM) 
DCQD 



</-" 



+ 1 



1 

l(d +!) + !]" 



1-body 
joint 1-body 
MUB 
POVM 
Stabilizer/Normalizer 



a 7i,: the Hilbert space of each experimental configuration 
c overall number of experimental configurations 



single-body 
single-body 
many -body 
many-body 
single- and two- body 



• • <g> X qm -r ZP™r I34D. Therefore, we have: 

N 2 -l 



£(p) = J2 Xmn(x q ^zp^®...®x q ^z p ^)p 1 ®...®p n (x^zp^^...®x q ^z p ^y 



m,n— 

d 2 -i 



^2 X( mi ... mT )(n 1 ...n r ){E mi PiE\ i )® . . . (E ms p s El s ). . . <E> (E mr p r El r ) 

mi,...,m r ,fii. , . . .,n r —Q 



d?-X 

= X] X(m 1 ...m r ){n 1 ...n r )(E m pEl)f , 

rci\...m T .n\. ..n r —0 

where we have introduced E nu = X qm ° Z Pm ° and Xmn = X(m 1 ,...,m r )(n 1 ....,n r )- I-e., m = (mi, . . . ,m s , . . . , m r ) and 
n = (ni, . . . , n s , . . . , n r ), and the index s represents a generic qudit. Let us first investigate the tensor product structure of the 
DCQD algorithm for characterization of the diagonal elements of the superoperator. We prepare the eigenstate of the stabilizer 
operators (Ef-E^)f and (E^E^ )f . For each qudit, the projection operators corresponding to outcomes ui k and ui k (where 
k,k' = 0,1,. ..,d-l), have the form P k = \ u~ lk {EfEf) 1 m&P k , = ± Ya~=o UJ ^ l ' k ' {E$Ef) v . The joint probability 

distribution of the commuting Hermitian operators P kl , , Pk 2 > Pk' 2 >■■■■> Pk r , Pk' r on tne output state £(p) is: 

- d 2 -l 

Tr[(P k P k ,)fS(p)] = X( mu ...,mr)(n lt ...,n r ) X 

mi,...,m r ,ni,. . .,n r — 

/d-id-i \ ® r 



Y,Y, u ~ lku ~ l ' k '^ E n{ E t)\E^) 1 ' E m {Ef) l {Ef,f p] 



\i=0 J'=0 / s 

By introducing EiE m = w lm E m Ei for each qudit and using the relation [{EfEf) l (E$Ef, ) l 'p] a = rt s we obtain: 

d 2 -i /d-id-i \® r 

Tr[(P k P k ,)fE( P )] = _ X(m u ..,m r ) { n u ..,n v )X ^ ^uMV^^] 

^ ' mi,...,m r ,ni,...,n r =G \i=0 i'=0 / a 

Using the QEC condition for nondegenerate codes, Tr[E^E m p] s = (5 mn ) s , and also using the discrete Fourier transform 
identities Y^iZo u^ m ~ k ^ 1 = d6i m , k and Xw'=o ^ lm ^ k ^ = dS? tk i for each qudit, we get: 

d 2 -i 

Tr[(P k P k ,)f£(p)] = X(m 1 ,..., mr )(nx,...,n.)(^,fc^,fc'^r l )f 
mi,...,m n tti,.. . ,n r — 
— X(moi,...,m 0r )(moi ; ...,mi) r ), 

where for each qudit, the index m 1S defined through the relations i mo = k and i' m = k', etc. I.e., E mo is the unique error 
operator that anticommutes with the stabilizer operators of each qudit with a fixed pair of eigenvalues u> k and u> k corresponding 
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to experimental outcomes k and k'. Since Pk and P).i operator have d eigenvalues, we have d 2 possible outcomes for each 
qudit, which overall yields (d 2 ) r equations that can be used to characterize all the diagonal elements of the superoperator with 
a single ensemble measurement and (2d) r detectors. Note that in the above ensemble measurement we can obtain log 2 d ir bits 
of classical information, which is optimal according to the Holevo bound for an 2r-qudit system of dimension d 2 . Similarly, the 
off-diagonal elements of superoperators can be identified by a tensor product of the operations in the DCQD algorithm for each 
individual qudit, see Ref. ||26l]. A comparison of the required physical resources for n qudits is given in Table U 

For a d-dimensional quantum system where d is neither a prime nor a power of a prime, we can always imagine another d'- 
dimensional quantum system such that d' is prime, and embed the principal qudit as a subspace into that system. For example, 
the energy levels of a six-level quantum system can be always regarded as the first six energy levels of a virtual seven-level 
quantum system, such that the matrix elements for coupling to the seventh level are practically zero. Then, by considering the 
algorithm for characterization of the virtual seven-level system, we can perform only the measurements required to characterize 
superoperator elements associated with the first six energy levels. 
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